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I. INTRODUCTION 

The increasing bulk of data that have been accumu- 
lated in the last few years have paved the way to the 
emergence of a new standard cosmological model usu- 
ally referred to as the concordance model. The Hubble 
diagram of Type la Supernovae (hereafter SNela), mea- 
sured by both the Supernova Cosmology Project [lj and 
the High-z Team p| up to redshift z ~ 1, first indi- 
cated that the universe is undergoing a phase of acceler- 
ated expansion. On the other hand, balloon born exper- 
iments, such as BOOMERanG and MAXIMA ||, de- 
termined the location of the first and second peak in the 
anisotropy spectrum of cosmic microwave background ra- 
diation (CMBR) strongly indicating that the geometry 
of the universe is spatially flat. If combined with con- 
straints coming from galaxy clusters on the matter den- 
sity parameter Omj these data indicate that the universe 
is dominated by a non-clustered fluid with negative pres- 
sure generically dubbed dark energy which is able, even- 
tually, both to close the universe and drive its accelerated 
expansion. This picture has been further strengthened by 
the more precise measurements of the CMBR spectrum, 
due to the WMAP experiment Q, and by the extension 
of the SNela Hubble diagram to z > 1 ■ 

After the discovery of this scenario, an overwhelm- 
ing flood of papers, presenting a great variety of mod- 
els which try to explain this phenomenon has appeared; 
in any case, the simplest explanation is claiming for the 
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well known cosmological constant A @. Although the 
best fit to most of the available astrophysical data [6j , the 
ACDM model failed in explaining why the inferred value 
of A is so tiny (120 orders of magnitude lower) compared 
to the typical vacuum energy values predicted by parti- 
cle physics and why its energy density is comparable to 
the matter density right now (the so called coincidence 
problem). As a tentative solution, many authors have 
replaced the cosmological constant with a scalar field 
rolling down its potential and giving rise to the model 
now referred to as quintessence |9(. Even if successful in 
fitting the data, the quintessence approach to dark energy 
is still plagued by the coincidence problem since the dark 
energy and matter densities evolve differently and reach 
comparable values for a very limited portion of the uni- 
verse evolution that is taken to coincide with the present 
era. Moreover, it is not clear where this scalar field origi- 
nates from thus leaving a great uncertainty on the choice 
of the scalar field potential (see [nj for comprehensive 
reviews) . 

The subtle and elusive nature of dark energy has led 
many authors to look for completely different scenarios 
able to give a quintessential behavior without the need 
of exotic components. 

To this aim, it is worth stressing that the acceleration 
of the universe only claims for a negative pressure dom- 
inant component, but does not tell anything about the 
nature and the number of cosmic fluids filling the uni- 
verse. This consideration suggests that it could be possi- 
ble to explain the accelerated expansion by introducing 
a single cosmic fluid with an equation of state causing it 
to act like dark matter at high densities and dark energy 
at low densities. An attractive feature of these mod- 
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els, usually referred to as Unified Dark Energy (UDE) or 
Unified Dark Matter (UDM) models, is that such an ap- 
proach naturally solves, al least phenomenologically, the 
coincidence problem. Some interesting examples are the 
generalized Chaplygin gas Ol, the tachyon field 01 an d 
the condensate cosmology \1$ • 

A different class of UDE models has been recently pro- 
posed by some of us 14] where a single fluid is considered 
whose energy density scales with the redshift in such a 
way that the radiation dominated era, the matter dom- 
ination era and the accelerating phase can be naturally 
achieved. It is worth noting that such a model is ex- 
tremely versatile since it can be interpreted both in the 
framework of UDE models and as a two-fluid scenario 
with dark matter and scalar field dark energy. 

Actually, there is still a different way to face the prob- 
lem of cosmic acceleration. As stressed in Lue et al. , 
it is possible that the observed acceleration is not the 
manifestation of another ingredient in the cosmic pie, 
but rather the first signal of a breakdown of our under- 
standing of the laws of gravitation. From this point of 
view, it is thus tempting to modify the Friedmann equa- 
tions to see whether it is possible to fit the astrophysical 
data with a model comprising only the standard mat- 
ter. Interesting examples of this kind are the Cardassian 
expansion ^(| and the DGP gravity |l7j . 

Moving in this same framework, it is possible to find 
alternative schemes where a quintessential behavior is ob- 
tained by simply taking into account the effective contri- 
bution to cosm ology o f some usually neglected funda- 
mental physics |l8t ll9L l2Cj . For instance, a cosmolog- 
ical constant term may be recovered as a consequence 
of a non - vanishing torsion field thus leading to a model 
which is consistent with both SNela Hubble diagram and 
Sunyaev- Zel'dovich data coming from clusters of galax- 
ies |21|. SNela data could also be efficiently fitted in- 
cluding higher-order curvature invariants in the gravity 
Lagrangian It is worth noting that these alterna- 

tive schemes provide naturally a cosmological component 
with negative pressure whose origin is simply related to 
the geometry of the universe thus overcoming the prob- 
lems linked to the physical significance of the scalar field. 

It is evident, from this short overview, the high num- 
ber of cosmological models which are viable candidates to 
explain the observed accelerated expansion. This abun- 
dance of models is somehow a result of having only a lim- 
ited number of cosmological tests to discriminate among 
rival theories so that any new observable which could be 
added to the usual ones is welcome. To this aim, it is 
useful to remember that both the SNela Hubble diagram 
and the recently proposed angular size - redshift relation 
of compact radio sources [2^ are distance based methods 
to probe cosmological models. From this point of view, 
it is interesting to look for tests based on time-dependent 
observables. To this aim, one can take into account the 
lookback time to distant objects since this quantity can 
discriminate among different cosmological models. The 
lookback time is observationally estimated as the differ- 



ence between the present day age of the universe and the 
age of a given object at redshift z. Such an estimate is 
possible if the object is a galaxy observed in more than 
one photometric band since its color is determined by its 
age as a consequence of stellar evolution. It is thus pos- 
sible to get an estimate of the galaxy age by measuring 
its magnitude in different bands and then using stellar 
evolutionary codes to choose the model that reproduces 
the observed colors at best. A q uite similar approach was 
pursued by Lima & Alcaniz |24| who used the age (rather 
than the lookback time) of old high redshift galaxies to 
constrain the dark energy equation of state (see also pjjl. 
The same method was then applied also to braneworld 
models [2(j and the Chaplygin gas p7j . 

It is worth noting, however, that the estimate of the 
age of a single galaxy may be affected by systematic er- 
rors which are difficult to control. Actually, this problem 
can be overcome by considering a sample of galaxies be- 
longing to the same cluster. In this way, by averaging 
the estimates of all galaxies, one obtains an estimate of 
the cluster age and reduces the systematic errors. Such 
a method was first proposed by Dalai et al. and 
then used by Ferreras et al. [23 to test a class of mod- 
els where a scalar field is coupled with the matter term 
giving rise to a particular quintessence scheme. We im- 
prove here this analysis using a different cluster sample 
and testing three different approaches to the dark energy 
problem. Moreover, we add a further constraint to bet- 
ter test the dark energy models and assume that the age 
of the universe for each model is in agreement with re- 
cent estimates. Note that this is not equivalent to the 
lookback time as we will discuss below. 

The layout of the paper is the following. In Sect. II, 
we briefly present the three classes of cosmological mod- 
els considered giving also the main quantities which we 
need for lookback time test. The method is developed in 
Sect. Ill and the used data are presented in Sect. IV. A 
detailed discussion of the results is the subject of Sect. V. 
We summarize and conclude in Sect. VI. 



II. DARK ENERGY MODELS 

As overviewed in the Introduction, many rival theories 
have been advocated to fit the observed accelerated ex- 
pansion and to solve the puzzle of the nature of the dark 
energy. Actually, we may divide the different cosmologi- 
cal models in three wide classes. According to the models 
of the first class, the dark energy is new ingredient of the 
cosmic Hubble flow, the simplest case being the ACDM 
scenario and its generalization which we will refer to as 
QCDM model. This is in sharp contrast with the as- 
sumption of UDE models (the second class) where there 
is a single fluid described by an equation of state compre- 
hensive of all regimes of cosmic evolution such as the one 
proposed in |L4| which we will consider here referring to 
it as the parametric density model. Finally, according to 
the third class models, accelerated expansion is the first 
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evidence of a breakdown of the Einstein General Relativ- 
ity (and thus of the Friedmann equations) which has to 
be considered as a particular case of a more general the- 
ory of gravity. As an example of this kind of models, we 
will consider the /(i?) -gravity [H, [H HJ H2 . Far from 
being exhaustive, considering these three models allow to 
explore very different scenarios proposed to explain the 
observed cosmic acceleration. 

In the following we will sketch these three dark energy 
models and derive all the quantities we will need for the 
analysis developed in Sect. III. 

A. The QCDM model 



The SNela Hubble diagram, the large scale galaxy clus- 
tering and the CMBR anisotropy spectrum can all be 
fitted by the ACDM model with (O M>0 , £l Q ) ~ (0.3, 0.7) 
thus giving qo ~ —0.55, i.e. the universe turns out to be 
in an accelerated expansion phase. The simplicity of the 
model and its capability of fitting the most of the data are 
the reasons why the ACDM scenario is the leading can- 
didate to explain the dark energy cosmology. Nonethe- 
less, we will first consider its generalization, the QCDM 
model, before concentrating on the ACDM scenario. 



Cosmological constant A has now become a textbook 
candidate to drive the accelerated expansion of the spa- 
tially flat universe. Despite its conceptual problems, the 
ACDM model turns out to be the best fit to a com- 
bined analysis of completely different astrophysical data 
ranging from SNela to CMBR anisotropy spectrum and 
galaxy clustering [6J, l30l l3l| . As a simple generaliza- 
tion, one may consider the QCDM scenario in which the 
barotropic factor w = p/p takes a negative value with 
w = — 1 corresponding to the standard cosmological con- 
stant. How such a negative pressure fluid drives the cos- 
mic acceleration may be easily understood by looking at 
the Friedmann equations : 
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where the dot denotes differentiation with respect to 
cosmic time t, H is the Hubble parameter and we 
have assumed that the universe is spatially flat as sug- 
gested by the position of the first peak in the CMBR 
anisotropy spectrum 0. 0. From the continuity equa- 
tion, p + 3H(p + p) = 0, we get for the i-th fluid with 

Pi =WiPi\ 



B. The parametric density model 

In the framework of UDE models, it has been recently 
proposed 0] a phenomenological class of models by in- 
troducing a single fluid 1 with energy density : 
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with < a < (3, s and b (with s < b) two scaling fac- 
tors and A norm a normalization constant. For several 
applications, it is useful to rewrite the energy density as 



a function of the redshift z. 
Eq.©, we get : 



Replacing a = (1 + z) in 



1 + z 
l + z s 



0-a 



1 + Z 
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(7) 



having defined z s = 1/s — 1 and Zb = 1/b— 1. It is easy 
to show that p cx aT® for a « s, p oc a~ a for s << 
a << b and p cx const, for a >> b. By setting (a, (3) = 
(3, 4) the energy density smoothly interpolates from a 
radiation dominated phase to a matter dominated period 
finally approaching a de Sitter state. The normalization 
constant A norm may be estimated by inserting Eq.© 
into Eq. Q and evaluating the result today. This gives : 



^, a~ 3(1+tUl) =^,o(l + ^) 3(1+lUl) , (3) 



where z = l/a — 1 is the redshift, fl; = pi/p cr it is the den- 
sity parameter for the i-th fluid in terms of the critical 
density p cr it — 3Hq/8itG and, hereafter, we label all the 
quantities evaluated today with a subscript 0. Inserting 
this result into Eq. , we get : 



Pcrit.O 



(1 + s)/ 3 -" (l + b a ) ' 



(8) 



The continuity equation may be recast in a form that 
allows us to compute the pressure for our fluid and then 
the barotropic factor w — p/p obtaining 0] : 



H{z) = H y/n M ,a(l + z f + "q,o(1 + z)W+") . (4) 

Using Eqs.QJ, J5J and the definition of the deceleration 
parameter q = —aa/a 2 , one finds : 



It is worth stressing that this model may be interpreted not only 
as comprising a single fluid with an exotic equation of state, but 
also as made of dark matter and scalar field dark energy or in 
the framework of modified Friedmann equations. Here, we prefer 
the UDE interpretation even if the results do not depend on this 
assumption. 
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[(a - 3)o + f/3 - 3)s]6 tt - [3 (a + s) + (a - [3)s]a a 

w = 

3(a + s)(a a + b a ) 

(9) 

which shows that the barotropic factor strongly depends 
on the scale factor (and hence on the redshift). Combin- 
ing the Friedmann equations, we get for the deceleration 
parameter q = (1 + 3u>)/2 which, in our case, gives : 



_ [(a - 2)a + {f3- 2)s}b a - [2(o + s) + {a - (3)s]a a 
q ~ 2(a + s){a a + b a ) 

inserting a = 1 gives the present day value as : 



(10) 



(y - l)q + z s [a y - 2(1 + y)} + ((3- 4)(1 + y) 
2(2 + z s )(l + y) 



(11) 

with y = (1+Zb) a . Some straightforward considerations 
allows to derive the following constraints on qo 0] : 
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2(2 + z s ) 



(12) 



It is convenient to solve Ea. (|ll|l with respect to Zb in 
order to express this one as a function of qo and z s . It 
is : 



a(l + z s ) + (3 - {2 + z s )(2g a + 2) 
a- (3+(2 + z s )(2q + 2) 



l/c 



1 • (13) 



The parametric density model is fully characterized by 
five parameters which are chosen to be the two asymp- 
totic slopes (a, (3), the present day values of the decel- 
eration parameter and of the Hubble constant (qo,Ho) 
and the scaling redshift z s . As in [l4|, we will set 
(a, (3) = (3,4) and z s = 3454 so that (qo,H ) will be 
the two parameters to be constrained by the data. 



C. Curvature quintessence 



where the energy density and the pressure due to the 
higher - order curvature invariants are [l8L |20| : 



1 



f'{R) 



f(R) - Rf'(R)] - 3HRf"(R) 



(14) 



/7 i^ {{2HR + R)r{R) + R 2 f"{R) 



-~[f{R)-Rf{R)\ 



(15) 



This approach is particularly useful since it allows to in- 
terpret the non-Einstein part of gravitational interaction 
as a new fluid with energy density and pressure given by 
(|14fl and (|15|) respectively. The barotropic factor for such 
a fluid turns out to be : 



«v 



-1 



f"(R)R+ f"(R)R-Hf"(R) 



[f(R)-Rf'(R)}/2-3Hf"{R) 



(16) 



A leading role is played by the form chosen for f(R). Fol- 
lowing [lSf. we set f(R) = foR n with /o constant. Unfor- 
tunately, it is not possible to analitically solve the Fried- 
mann equations for the model with both matter and cur- 
vature contributions. However, if we neglect the matter 
term 2 (i.e. we set pm = 0), it is possible to find power law 
solutions of the Friedmann equations as a(t) — (t/to) a 
with US : 



2n 2 -371+1 



(17) 



This gives an entire class 3 of cosmological models with 
constant deceleration parameter : 



q(t) = qo = 



1 - a 



2n 2 - 2n - 1 
2n 2 - 3n + 1 



(18) 



In order to get an accelerated expansion (a > and 
qo < 0), the parameter n has to satisfy the following 
constraint : 



There is no a priori reason to restrict the gravity La- 
grangian to a linear function of Ricci scalar R as in the 
usual Einstein general relativity. Actually, higher order 
curvature invariants are introduced in quantum gravity 
theories so that it is worth considering the effect of such 
generalizations on the late evolution of the universe. 

Replacing the Ricci scalar R with a (up to now) generic 
function f(R) of Ricci scalar in the gravity lagrangian, 
the resulting field equations may still be recast in the 
Friedmann - like form provided that the total energy den- 
sity and the total pressure are written as : 



Ptot = Pm + P CUTV 



Ptot = Pm + Pc 



n € 



-00, ■ 



1 - V3 



1 + V3 



(19) 



The model we obtain by this ansatz should be considered as a 
toy model that is able to give some useful insights on the true 
model. However, it is worth noting that the matter contribution 
is indeed much smaller than the curvature one and it should 
be indeed negligible if only the baryonic matter is considered 
provided that the theory is able to erase the need for dark matter. 
Actually, there is another class of models characterized by a = 
and n = 0, 1/2, 1, but it is ruled out by the non vanishing value 
of the Hubble constant. 
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We will refer to models with n in the first (second) range 
as CurvDown (CurvUp) models respectively and use the 
method described in the following section to constrain 
the two parameters (n, Hq) which completely assign the 
model. 



III. THE METHOD 



j-ob. 



(Zi) = t L (z F ) - ti{z) 



t L (z F )] 



(22) 



where t% bs 



is the today estimated age of the universe and 
a delay factor can be defined as : 



Most of the tests recently used to constrain cosmolog- 
ical parameters (such as the SNela Hubble diagram and 
the angular size - redshift) are essentially distance - based 
methods. The proposal of Dalai et al. [HJ to use the 
lookback time to high redshift objects is thus particu- 
larly interesting since it relies on a completely different 
observable. The lookback time is defined as the differ- 
ence between the present day age of the universe and its 
age at redshift z and may be computed as : 



*zO,p) = t H 



dz' 



(l + z')E(z',p) 



(20) 



where t H = 1/H = 9.78/1" 1 Gyr is the Hub- 
ble time (with h the Hubble constant in units of 
100 km s _1 Mpc" 1 ), E(z,p) = H(z)/H is the dimen- 
sionless Hubble parameter and we have denoted with {p} 
a set of parameters characterizing a given cosmological 
model. It is worth noting that, by definition, the look- 
back time is not sensible to the present day age of the 
universe to so that it is (at least in principle) possible 
that a model fits well the data on the lookback time, but 
nonetheless it predicts a completely wrong value for t^. 
This latter parameter can be evaluated from Eg. 12011 by 
simply changing the upper integration limit from z to in- 
finity. This shows that it is indeed a different quantity 
since it depends on the full evolution of the universe and 
not only on how the universe evolves from the redshift 
z to now. That is why we will explicitly introduce this 
quantity as a further constraint. 

Let us now discuss how we use the lookback time and 
the age of the universe to test a given cosmological model. 
To this end, let us consider an object i at redshift z and 
denote by ti{z) its age defined as the difference between 
the age of the universe when the object was born, i.e. at 
the formation redshift zp, and the one at 2. It is : 



U(z) 



dz' 



dz' 



(l + z')E(z',p) J ZF (l + z')E(z',p) 
dz' 



(l + z')E(z',p) 



(21) 



where, in the last row, we have used the definition (|20|l of 
the lookback time. Suppose now we have N objects and 
we have been able to estimate the age of the object 
at redshift zi for i = 1,2, ...,N. Using the previous 
relation, we can estimate the lookback time f^ s (zi) as : 



df = t° bs -t L (z F ) 



(23) 



The delay factor is introduced to take into account our 
ignorance of the formation redshift zp of the object. Ac- 
tually, what can be measured is the age U(z) of the object 
at redshift z. To estimate zp, one should use Ea. l|21|) 
assuming a background cosmological model. Since our 
aim is to determine what is the background cosmolog- 
ical model, it is clear that we cannot infer zp from the 
measured age so that this quantity is completely undeter- 
mined. It is worth stressing that, in principle, df should 
be different for each object in the sample unless there is a 
theoretical reason to assume the same redshift at the for- 
mation of all the objects. If this is indeed the case (as we 
will assume later), then it is computationally convenient 
to consider df rather than zp as the unknown parameter 
to be determined from the data. We may then define a 
likelihood function as : 



Cu(p, df) oc exp [-xtt(P, df)/2] 



(24) 



with : 



Xit 



N-N p + 1 



4-theor ( ~\ j-o 
z \P) ~ z 



obs 



N 

E 



tf e ° r ( Zi ,p) 



4-obs 



(Zi) 



(25) 



where N p is the number of parameters of the model, at 
is the uncertainty on ig hs , Oj the one on t L bs (zi) and the 
superscript theor denotes the predicted values of a given 
quantity. Note that the delay factor enters the defintion 
of xft smcc it determines t c j bs (zi) from ti(z) in virtue 
of Ea. (|22(l . but the theoretical lookback time does not 
depend on df. 

In principle, such a method should work efficiently to 
discriminate among the various dark energy models. Ac- 
tually, this is not exactly the case due to the paucity of 
the available data which leads to large uncertainties on 
the estimated parameters. In order to partially alleviate 
this problem, it is convenient to add further constraints 
on the model by using a Gaussian prior 4 on the Hubble 
constant, i.e. redefining the likelihood function as : 



4 The need for a prior to reduce the parameter uncertainties has 
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£(p) cx £; t (p)exp 



1 fh-h obs 

2 l _ V h 



cx exp[-x 2 (p)/2] 
(26) 



where we have absorbed df in the set of parameters p 
and have defined : 



X =Xit + 



h-h 



ob.h 



(27) 



with h obs the estimated value of h and ah its uncertainty. 
We use the HST Key project results [33 setting (h, ah) — 
(0.72, 0.08). Note that this estimate is independent of the 
cosmological model since it has been obtained from local 
distance ladder methods. 

The best fit model parameters p may be obtained by 
maximizing £(p) which is equivalent to minimize the x 2 
defined in Ea. (|27|l . It is worth stressing that such a func- 
tion should not be considered as a statistical \ 2 i n the 
sense that it is not forced to be of order 1 for the best 
fit model to consider a fit as successful. Actually, such 
an interpretation is not possible since the errors on the 
measured quantities (both £j and to) are not Gaussian 
distributed and, moreover, there are uncontrolled sys- 
tematic uncertainties that may also dominate the error 
budget. Nonetheless, a qualitative comparison among 
different models may be obtained by comparing the val- 
ues of this pseudo \ 2 even if this should not be considered 
as a definitive evidence against a given model. 

Having more than one parameter, we obtain the best 
fit value of each single parameter pi as the value which 
maximizes the marginalized likelihood for that parameter 
defined as : 



C Pi cx J dpi . . . J dpi-i J dp l+x . . .J dp n £(p) . (28) 

After having normalized the marginalized likelihood to 1 
at maximum, we compute the la and 2a confidence limits 
(CL) on that parameter by solving C Pi = cxp (—0.5) and 
C Pi = exp (—2) respectively. 



in this sense since they can be detected up to high red- 
shift and their redshift at formation epoch 5 is almost 
the same for all the clusters. Furthermore, it is relatively 
easy to estimate their age from photometric data only. 

To this end, the color of their component galaxies, in 
particular the reddest ones, is all what is needed. Ac- 
tually, the stellar populations of the reddest galaxies be- 
come redder and redder as they evolve. It is just a mat- 
ter, then, to assume a stellar population synthesis model, 
and to look at how old the latest episode of star forma- 
tion should be happened in the galaxy past to produce 
colors as red as the observed ones. This is what we will 
refer to as color age. The main limitation of the method 
relies in the stellar population synthesis model, and on a 
few (unknown) ingredients (among which the metallicity 
and the star formation rate law). 

The choice of the evolutionary model is a key step in 
the estimate of the color age and the main source of un- 
certainty |33|. An alternative and more robust route to 
cluster age is to consider the color scatter (see [34| for an 
early application of this approach). The argument, qual- 
itatively, goes in this way : if galaxies have an extreme 
similarity in their color and nothing is conspiring to make 
the color scatter surreptitiously small, then the latest 
episode of star formation should happened in the galaxy 
far past, otherwise the observed color scatter would be 
larger. Quantitatively, the scatter in color should thus 
be equal to the derivative of color with time multiplied 
the scatter of star formation times. The first quantity 
may be predicted using population synthesis models and 
turns out to be almost the same for all the evolutionary 
models thus significantly reducing the systematic uncer- 
tainty. We will refer to the age estimated by this method 
as scatter age. 

The dataset we need to apply the method described 
in the previous section may now be obtained using the 
following procedure. First, for a given redshift Zj, we 
collect the colors of the redddest galaxies in a cluster at 
that redshift and then use one of the two methods out- 
lined above to determine the color or the scatter age of 
the cluster. If more than one cluster is available at that 
redshift, we average the results from different clusters in 
order to reduce systematic errors. Having thus obtained 
ti(zi), we then use Ea. (|22|) to estimate the value of the 



IV. THE DATA 



In order to apply the method outlined above, we need 
a set of distant objects whose age can be somehow esti- 
mated. Clusters of galaxies seem to be ideal candidates 



been also advocated in previous works using the age of old high 
redshift galaxies as a cosmological test. For instance, in Ref. |24| 
a strong prior on Qm is introduced to better constrain the dark 
energy equation of state. It is likely that extending the dataset 
to higher redshifts and reducing the uncertainties on the age 
estimate will allow to avoid resorting to priors. 



5 It is worth stressing that, in literature, the cluster formation 
redshift is defined as the redshift at which the last episode of 
star formation happened. In this sense, we should modify our 
definition of df by adding a constant term which takes care of 
how long is the star formation process and what is the time 
elapsed from the beginning of the universe to the birth of the 
first cluster of galaxies. For this reason, it is still possible to 
consider the delay factor to be the same for all clusters, but it 
is not possible to infer zp from the fitted value of df because 
we do not know the detail of the star formation history. We 
stress that this approach is particular useful since it allows us 
to overcome the problem to consider lower limits of the universe 
age at z rather than the actual values. 
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Color age 


Scatter age 


z 


N 


Age (Gyr) 


Ref 


z 


N 


Age (Gyr) 


Ref 


0.60 


1 


4.53 


[35] 


0.10 


55 


10.65 


f38] 


0.70 


3 


3.93 


[36] 


0.25 


103 


8.89 


[38J 


0.80 


2 


3.41 


f35] 


1.27 


1 


1.60 


f39] 



TABLE I: Main properties of the cluster sample used for the 
analysis. The data in the left part of the Table refers to 
clusters whose age has been estimated from the color of the 
reddest galaxies (color age) , while that of clusters in the right 
part has been obtained by the color scatter (scatter age). For 
each data point, we give the redshift z, the number TV of 
clusters used, the age estimate and the relevant reference. 



lookback time at that redshift. Actually, what we mea- 
sure is t^ 3 (zi) + df that is the quantity that enters the 
definition l|25l) of xft an d then the likelihood function. 

To estimate the color age, following [3^, |3(| , we have 
chosen, among the various available stellar population 
synthesis models, the Kodama and Arimoto one |37| . 
which, unlike other models, allows a chemical evolution 
neglected elsewhere. This gives us three points on the 
diagram z vs. i£ obtained by applying the method to a 
set of six clusters at three different redshifts as detailed in 
the left part of Table 1 . Using a large sample of low red- 
shift SDSS clusters imaged, one of us [11] has evaluated 
the scatter age for clusters age at z — 0.10 and z — 0.25, 
while Blakeslee et al. [3£j applied the same method to 
a single, high redshift (z = 1.27) cluster. Collecting the 
data using both the color age and the scatter age, we end 
up with a sample of ~ 160 clusters at six redshifts (listed 
in Table 1) which probe the redshift range (0.10,1.27). 
This nicely overlaps the one probed by SNela Hubble di- 
agram so that a comparison among our results and those 
from SNela is possible. We assume a a = 1 Gyr as uncer- 
tainty on the cluster age, no matter what is the method 
used to get that estimate. Note that this is a very con- 
servative choice. Actually, if the error on the age were 
so large, the color - magitude relation for reddest cluster 
galaxies should have a large scatter that is not observed. 
We have, however, chosen such a large error to take quali- 
tatively into account the systematic uncertainties related 
to the choice of the evolutionary model. 

Finally, we need an estimate of tg bs to apply the 
method. Following Rebolo et al. [^U, we choose 
(tQ bs ,at) — (14.4,1.4) Gyr as obtained by a combined 
analysis of the WMAP and VSA data on the CMBR 
anisotropy spectrum and SDSS galaxy clustering. Ac- 
tually, this estimate is model dependent since Rebolo et 
al. |31j implicitly assumes that the ACDM model is the 
correct one. However, this value is in perfect agreement 
with ig hs = 12.6]_2;4 Gyr determined from globular clus- 
ters ag e l40l and ig bs > 12.5±3.5 Gyr from radioisotopes 
studies |4lj . For this reason, we are confident that no sys- 
tematic error is induced on our method using the Rebolo 
et al. estimate for tg &s even when testing cosmological 



models other than the ACDM one. 



V. RESULTS 

We have applied the method outlined above to the dark 
energy models described in Sect. II in order to constrain 
their parameters and see if they are viable candidates 
to explain cosmic acceleration. To this aim, let us note 
that each one of the presented models is fully described 
by only few parameters which are 6 (^M,h, w) for the 
QCDM model, (qo,h) for the parametric density model 
and (n, h) for the curvature quintessence. For all three 
models, there is still another parameter entering the fit, 
that is the delay factor df, which we will marginalize over 
since it is not interesting for our aims. 

Let us first consider the QCDM model. The main re- 
sults are plotted in Figs.^ and |21 In the first plot, we 
compare the estimated clusters age with the quantity : 

t(z) = t L (z) + df (29) 

using the best fit values for the model parameters and 
the delay factor which turn out to be : 

{Q M , h, w) = (0.25,0.70,-0.81) , df = 4.5 Gyr 

giving x 2 — 0.04. The \ 2 value for the best fit parameters 
(both for the QCDM model and the other dark energy 
models considered) turns out to be quite small suggesting 
that we have seriously overestimated the errors. This is 
not surprising given the (somewhat arbitrary) way we 
have fixed the uncertainties on the estimated age of the 
clusters. That this is likely to be the case is also suggested 
by a qualitative argument. One could rescale the errors 
on t°£ s {zi) in such a way that x 2 = 1 f° r the best fit 
model. Since for the best fit QCDM model, x 2 — Xiv 
this leads to multiply by almost 1/5 the uncertainties on 
t°^ s {zi). If the error on to were negligible, this means that 
we should reduce the uncertainty on the cluster age from 
1 Gyr to 0.2 Gyr that is indeed a more realistic value. 
The presence of an error on £g bs slightly complicates this 
argument, but does not change the main conclusion. We 
are thus confident that the very low value of the x 2 we g e t 
for the best fit model is only due to overestimating the 
uncertainties on the clusters ages. However, we do not 
perform any rescaling of the uncertainties since, to this 
end, we should select a priori a model as the confidence 
one which is contrary to the philosophy of the paper. It 
is worth stressing, however, that such rescaling does not 
affect anyway the main results. 



6 We drop the subscript from !1m,o since it does not give rise to 
any confusion here. Also, we use the dimensionless parameter h 
instead of the Hubble constant Hq . 
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FIG. 1: Comparison among predicted and observed values of 
t = t L (z) + df for the best fit QCDM model. 
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FIG. 2: The la and 2a confidence regions in the plane 
(Qm,w) for the QCDM model after marginalizing over the 
Hubble constant h and the delay factor df. 



Fig-El shows the la and 2a confidence levels in the 
(Om, w) plane obtained by marginalizing over the Hubble 
constant and the delay factor. Two interesting consider- 
ations may be drawn from that plot. First, we note that 
also phantom models (i.e. models with w < — 1 violat- 
ing the weak energy condition) are allowed by the data. 
This is in agreement with recent results coming from fit- 
ting the QCDM model to the SNela Hubble diagram and 
the CMBR anisotropy spectrum 42] . Unfortunately, a 
direct comparison is not possible since the marginalized 
likelihood is too flat to get any constraints on w so that 
all the values in the range tested (—2 < w < 1/3) are 
well within the la CL. Although essentially due to the 
paucity of the data, this result is also a consequence of 



T 
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FIG. 3: Comparison among predicted and observed values of 
t = th{z) + df for the best fit ACDM model. 
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FIG. 4: The la and 2a confidence regions in the plane (Qm, h) 
for the ACDM model. 

having not used any prior on m as it is usually done in 
other analyses. By using the procedure described at the 
end of Sect. Ill, we obtain the following estimates for the 
other QCDM parameters : 

9, M G (0.13,0.39) , h e (0.63, 0.77) (la CL) , 

fl M £ (0.01,0.62) , he (0.56,0.84) (2a CL) . 

Given that we are able to give only weak constraints on 
the QCDM model, from now on we will only dedicate our 
attention to the case w = — 1, i.e. to the ACDM model 
and do not discuss anymore the results for the QCDM 
model. The best fit parameters for the cosmological con- 
stant model turn out to be : 



0.2 0.4 0.6 0.8 1 1.2 1.4 



FIG. 5: Comparison among predicted and observed values of 
t = tr,(z) + d/ for the best fit parametric density model. 



(fl M ,h) = (0.22,0.71) , d/ = 4.05Gyr ( X 2 0.09) 

that gives rise to the curve t(z) shown in Fig. [31 while 
Fig. ^reports the confidence regions in the (f^M, h) plane 
after marginalizing over the delay factor. From the 
marginalized likelihood functions, we get : 

CL M e (0.10,0.35) , h G (0.63, 0.78) (la CL) , 

Q M € (0.06,0.59) , ft € (0.56,0.85) (2a CL) . 

The ACDM model has been widely tested against a 
large set of different astrophysical data. This offers us the 
possibility to cross check both the model and the valid- 
ity of method. It is instructive, in this sense, to compare 
our results with those coming from the fit to the SNela 
Hubble diagram. For instance, Barris et al. Q used a 
set of 120 SNela up to z = 1.03 finding Qm = 0.33 as 
best fit value with a large uncertainty (not quoted ex- 
plicitly, but easy to see in their Fig. 12) in good agree- 
ment with our result. A more recent result has been 
obtained by Riess et al. Q which have used a SNela 
Hubble diagram extending up to z — 1.55 and have found 
fi&r = 0.291q q3 still in agreement with our result. The 
ACDM model has also been tested by means of the an- 
gular size - redshift relation. Using a catalog of ultracom- 
pact radio sources and taking carefully into account sys- 
tematic uncertainties and selection effects, Jackson [4(j 
has found n M = 0.24t°;°? in almost perfect agreement 
with our estimate. 

Neither Barris et al. Q nor Riess et al. Q quote 
a best fit value for h since this parameter is infinitely 
degenerate with the supernovae absolute magnitude M 
when fitting to the SNela Hubble diagram. Nonethe- 
less, SNela may still be used to determine h by fitting 
the linear Hubble law to the low rcdshidt (z < 0.1) 
SNela. Using this method, Daly & Djorgovski [43| have 
found h = 0.664±0.08 in good agreement with our re- 
sult. Moreover, it is worth noting that our estimate for h 



. 65 
. 6 



-0.9 -0.8 -0.7 -0.6 -0.5 -0.4 -0.3 

qo 

FIG. 6: The la and 2<r confidence regions in the plane (go, h) 
for the parametric density model. 

turns out to be in agreement with estimates coming from 
different methods such as various local standard candles 
[32|. the Sunyaev- Zel'dovich effect in galaxy clusters 
and time delays in multiply imaged quasars 0] . Finally, 
let us quote the results obtained by Tegmark et al. [3(j 
which have performed a combined fit of the ACDM model 
to both the WMAP data on the CMBR anisotropy spec- 
trum and the galaxy power spectrum measured by more 
than 200,000 galaxies surveyed by the SDSS collabora- 
tion. They find fi M = 0.30±0.04 and h = 0.70±° o ° o i in 
very good agreement with our results. 

Perhaps, the most interesting result of testing the 
ACDM model with our method is not the success of this 
model (since it has yet been shown by a lot of evidences), 
but the substantial agreement among our estimates of 
the parameters (fijw, h) and those coming from different 
kinds of data. This is quite encouraging since it is an 
important successful cross check and makes us confident 
about the results one could obtain by applying it to other 
cosmological models. 

Let us now examine the results obtained for the para- 
metric density model, pictorially shown in Figs .[S] and [5] 
The best fit model is obtained for the following values of 
the parameters (qo, h) and of the delay factor : 

(q ,h) = (-0.68,0.71) , d/ = 4.20Gyr ( X 2 =s 0.07) . 
Marginalizing over df, we get : 

q 6 (-0.81, -0.47) , h S (0.64, 0.78) (la CL) , 
qo e (-0.89,-0.24) , h e (0.58,0.85) (2a CL) . 
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FIG. 7: The la and 2a confidence regions in the plane (n, h) 
for curvature quintessence in the CurvUp regime. 

Note that the 2a CL on the go parameter has been trun- 
cated at the upper end since it formally extends to values 
higher than the physically acceptable one. 

In 0|, the parameters of this model have been con- 
strained by using both the SNela Hubble diagram and the 
angular size - redshift relation. In particular, fitting the 
model to the SNela Hubble diagram gives h = 0.64+^1, 
while the physically acceptable range for q$ turns out to 
be in agreement with the data for go = —0.42 as best 
fit value. The present day deceleration parameter qp is 
better constrained using the data listed in Jackson jig 
to perform the angular size - redshift test thus obtaining 
go = — 0.64^q'}2 UM ■ Both these results are in very good 
agreement with our estimates so that we conclude that 
the parametric density model is a viable candidate alter- 
native to the ACDM model. 

Finally, let us discuss the results for the curvature 
quintessence. Figs.[7| and [S] report the confidence re- 
gions in the plane (n, h) for the CurvUp and CurvDown 
regime respectively after marginalizing over the delay fac- 
tor. The first striking feature is that the contour plots 
are not closed so that the marginalized likelihood func- 
tion gives only an upper (lower) limit to the parameter n 
in the CurvUp (CurvDown) regime. Formally, we get the 
following estimates for the best fit values in the CurvUp 
and CurvDown regime respectively : 

(n, h) = (1.367, 0.71) , df = 4.80 Gyr ( X 2 ^ 0.23) , 

( n , h) = (-0.367, 0.74) , df = 4.80 Gyr { x 2 ^ 0.21) , 

but the best fit value for n actually lies outside the inves- 
tigated range for this parameter. Being the confidence 



0.85 



0.8 



. 75 
. 7 



0. 65 




-0.6 -0.55 -0.5 -0.45 -0.4 



n 

FIG. 8: The la and 2a confidence regions in the plane (n, h) 
for curvature quintessence in the CurvDown regime. 

regions open, it is meaningless to give constraints on h, 
but nonetheless it is possible to infer the following limits : 

n< 1.402 at la CL , n < 1.424 at 2a CL , 

n > -0.508 at la CL , n > -0.606 at 2a CL , 

where the first row refers to the CurvUp regime and the 
second to the CurvDown one. These limits do not contra- 
dict the r ang es determined by fitting to the SNela Hubble 
diagram [23], but seems to be quite non-realistic. Actu- 
ally, the deceleration parameter corresponding to these 
values of n is quite small (go ~ —0.01) thus contradicting 
the evidence in favor of an accelerating universe. More- 
over, the results in |22j | have been obtained by using an 
old sample of SNela, including some SNela that have now 
been discarded from the Gold set in Riess et al. 0- On 
the other hand, it is worth noting that fitting a power 
law scale factor to the angular size - redshift relation for 
compact radio sources gives a ~ 1 |47j which, by using 
Ea.l|17jl. translates in an estimate for n that is in good 
agreement with our result. 

VI. CONCLUSIONS 

The impressive amount of data indicating a spatially 
flat universe in accelerated expansion has posed the prob- 
lem of dark energy and stimulated the search for cosmo- 
logical models able to explain such unexpected behav- 
ior. Many rival theories have been proposed to solve 
the puzzle of the nature of dark energy ranging from a 
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rolling scalar field to a unified picture where a single ex- 
otic fluid accounts for the whole dark sector (dark mat- 
ter and dark energy). Moreover, modifications of the 
gravity Lagrangian has also been advocated. Although 
deeply different in their underlying physics, all these sce- 
narios share the common feature of well reproducing the 
available astrophysical data. It is worth stressing, how- 
ever, that the most widely used cosmological tests and 
in particular the SNela Hubble diagram and the angular 
size - redhisft relation) are essentially based on distance 
measurements to high redshift objects and are thus af- 
fected by similar systematic errors. It is hence partic- 
ularly interesting to look for methods which are related 
to the estimates of different quantities. Being affected 
by other kinds of observational problems, such methods 
could be considered as cross checks for the results ob- 
tained by the usual tests and they should represent com- 
plementary probes for cosmological models. 

The technique we have devised here is a first step in 
this direction. We have used the present age of the uni- 
verse and the lookback time to galaxy clusters to build 
a sort of time diagram of the universe in order to recon- 
struct its age evolution. Relying on stellar evolutionary 
codes, the estimate of the lookback time is related to a 
different astrophysics than the distance based methods 
and it is thus free from any problem connected with the 
evolution of standard candles (such as the SNela absolute 
magnitude and the intrinsic linear size of radio sources). 
Actually, this technique could be affected by its own sys- 
tematics (such as, for instance, the degeneracy between 
age and metallicity), but these may be better controlled. 
Moreover, comparing the results thus obtained with those 
coming from distance based methods allows to strengthen 
the conclusions suggested by both techniques. 

Motivated by these considerations, we have first ap- 
plied our method to the ACDM model in order to esti- 
mate the present day values of the matter density pa- 
rameter fijvf and of the dimensionlcss Hubble constant h 
obtaining (at the la CL) : 



These values are in agreement with the previous esti- 
mates in literature. This can be considered an indepen- 
dent confirmation not only of the viability of the ACDM 
model, but also of our method. It is worth noting that 
the ACDM scenario receives further support from this 
test and henceforth the cosmological constant A still re- 
mains the best candidate to explain the dark energy puz- 
zle from an observational point of view. 

Nonetheless, the ACDM model is severely affected by 
conceptual problems so that it is worthwhile to look fo 
alternative approaches. This has stimulated a plenty of 
models where the cosmic acceleration is due to a dom- 
inant scalar field rolling down its potential. However, 
such a scheme, on the one hand, still does not solve the 
coincidence problem and, on the other hand, is plagued 
by the unidentified nature of the scalar field itself and the 



ignorance of its self interaction potential. These consid- 
erations have opened the way to different models that are 
able to give an accelerated expansion without the need of 
scalar fields. Two of these approaches have been tested 
in this paper. 

In |l4j . it has been proposed a phenomenological uni- 
fied model where a single fluid with a given energy density 
assigned by few parameters is able to fit both the SNela 
Hubble diagram and the angular size - redshift relation 
for ultracompact radio sources without the need of any 
cosmological constant or scalar field. This parametric 
density model has been tested here with lookback time 
method and the following estimates (at the la CL) for 
its characterizing parameters has been obtained : 

q a = -0.68±°;^ , h = 0.71±0.07 , 

in agreement with the results from the test performed in 
[l4j. It is worth stressing that our method turns out to 
be more efficient than the usual fit to the SNela Hubble 
diagram for the parametric density model. It is also in- 
teresting to note that both q and t predicted by this 
model are in almost perfect agreement with those com- 
puted for the ACDM model with the parameters (11m, h) 
as discussed above. This is not surprising since Eq.JfjJ) 
shows that, nowadays, the energy density of the para- 
metric density model is very similar to that of the ACDM 
model so that these two scenarios share most of the ob- 
servable properties referred to today quantities. How- 
ever, this does by no way mean that the two models are 
the same. Actually, it is the underlying philosophy that 
is completely different having now a single fluid rather 
than a cosmological constant dominating over the mat- 
ter term. Both models predict similar values ofthe today 
observed quantities simply because they are tied to re- 
produce the same data and not because they are two 
different manifestations of the same underlying physics. 

Another possible approach to the cosmic acceleration 
is to consider this feature as the first signal of the break- 
down of the Einstein General Relativity at some charac- 
teristic scale. In this picture, the universe is still dom- 
inated by standard matter, but the Friedmann equa- 
tions have to be modified as consequence of a different 
gravity Lagrangian. This philosophy inspired curvature 
quintessence scenarios where an effective dark energy 
is related to the properties of the function f(R) which 
replaces the Ricci scalar R in the gravity Lagrangian 
[Til IT^. |2(j. The coupling between matter and curva- 
ture for a theory with f(R) — foR n gives rise to fourth 
order nonlinear differential equations for the evolution of 
the scale factor that is not analytically solvable. This 
difficulty disappears if one considers a toy model with 
£Im = so that power law solutions, a ~ t a , are pos- 
sible provided that a is linked to n by Ea. ((T7|l . Such a 
model is particularly attractive from a theoretical view- 
point (since it allows to give a purely geometric inter- 
pretation of the dark energy) and has also been shown 
to fit well SNela Hubble diagram [23. We have shown 
here that, while it is able to pass successfully the loock- 
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back time test implemented in this paper, only weak con- 
straints can be imposed on the parameter n with the best 
fit value lying in a region corresponding to decelerating 
rather than accelerating models. We are thus tempted to 
conclude that this scenario could be rejected with some 
confidence. Actually, this conclusion is also suggested 
by the recent result by Riess et al. which, using the 
SNela Hubble diagram up to z = 1.55, have convincingly 
detected a change in the sign of the deceleration param- 
eter q. The power law solutions predict, on contrary, a 
constant q so that they are ruled out by the result in Riess 
et al. 7]. It is important to stress, however, that this is 
not a definitive exclusion of the curvature quintessence 
scenario. Actually, we have only considered power law 
solutions in absence of matter. It is indeed conceivable 
that including a matter term completely change the evo- 
lution of the scale factor possibly giving a past decelera- 
tion followed by the present acceleration thus leading the 
curvature quintessence model in agreement with what is 
suggested by the extended SNela Hubble diagram (see 
also |19| and references therein for a further discussion of 
curvature models). 

Having tested three different models, it is worth ask- 
ing what is the better one. Unfortunately, this is not 
possible on the basis of the test results only. Compar- 
ing the x 2 values, one could naively conclude that the 
parametric density model is the best one since it gives 
the lowest x 2 ■ However, while the difference among the 
X 2 values is significant between the parametric density 
model and the curvature quintessence scenario (almost 
an order of magnitude), it is too small (0.07 vs. 0.09) to 
conclude that the ACDM model is disfavored. We have 
thus to conclude that this test alone is not able to dis- 
criminate among these two dark energy candidates. On 
the other hand, Sandvik et al. 01 recently claimed that 
UDE models are not viable because the growth of den- 
sity perturbations will lead to matter power spectrum in 



disagreement with what is observed. This should be an 
evidence against the parametric density model. However, 
it is worth noting that Sandvik et al. explicitely consider 
the generalized Chaplygin gas model which is character- 
ized by a negative squared sound speed that seems to 
be the main cause of the anomalous growth of perturba- 
tions. For the parametric density model, the sound speed 
is always positive definite so that it is likely that the ar- 
gument of Sandvik et al. should be at least reconsidered. 

Finally, we would briefly comment on the possibility to 
ameliorate our method by increasing the maximum red- 
shift probed. To this aim, galaxy clusters do not appear 
as good candidates since it is quite difficult to detect a 
significative number of member galaxies up to redshifts 
larger than ~ 1.3. However, high redshift galaxies may 
be taken into account provided that they are detected 
in as more photometric bands as possible. This latter 
requirement is fundamental since it allows not only to 
better estimate the age of the galaxy, but also to infer 
constraints on its formation redshift zp that cannot be 
assumed to be the same for all the galaxies (as it has been 
possible for clusters). The i drop out technique allows to 
discover galaxies up to redshift z ~ 6 and thus let us 
hope to measure lookback time up to such high redshift. 
With its ability of both furnishing multicolor photome- 
try of high redshift galaxies (and thus better estimates 
of their color ages and formation redshift) and increasing 
the number of SNela, the GOODS survey [5(j seems to 
be the most promising source of cosmological constraints 
in the near future. 
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